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ATTORNEY DOCKET NO. 60.1564 

SPECTROSCOPIC PH MEASUREMENT USING OPTIMIZED MIXTURES 
OF REAGENTS TO EXTEND MEASUREMENT RANGE 

[0001] The present invention claims priority to co-pending United States Patent 
Application Serial No. 10/667,639 filed September 22, 2003, entitled "Determining 
Fluid Chemistry of Formation Fluid by Downhole Reagent Injection Spectral 
Analysis," which in turn claims priority to co-pending GB patent application 
no. 0227267.2 filed November 22, 2002, entitled "Apparatus and Method for 
Analysing Downhole Water Chemistry." Both of these patent applications are 
incorporated by reference herein in their entireties. 

FIELD OF THE INVENTION 
[0002] The present invention relates to pH measurement of aqueous samples using 
pH sensitive reagents (such as dyes) and, more particularly, to a reagent mixture that 
allows accurate pH measurement over an extended measurement range. 

BACKGROUND 

[0003] Spectroscopic (or colorimetric) method for pH measurement is a well- 
established technique, wherein a pH sensitive reagent (or dye) changes color based on 
the pH of the solution. The theory of spectroscopic pH measurement using a single 
reagent indicator is well known (see Bates, "Determination of pH: Theory and 
Practice" Chapter 6 (1964), incorporated by reference herein in its entirety) and can 
be characterized by the following equation: 

pH = pKa + log^ + log— (1) 
Ya ^ 

where Ka is the thermodynamic acid dissociation constant for the reagent, which is a 
function of temperature and pressure; A, B are concentrations of the acid and base 
forms of the reagent, respectively; y^, ye are activity coefficients of the acid and base 
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forms of the reagent, respectively, which are a function of temperature, pressure and 
ionic strength of the solution. 

[0004] The ratio of the base form to the acid form ifi/A) of the reagent indicator may 
be calculated from spectral measurements using the following equations: 

Cj=A + B (2) 
OD'" =ej'M + f^'/i5 (3) 
OD^^ =eiHA + eiHB (4) 



/ f \ 



(6) 



where OD^* is the optical density measured at wavelength A/; / is path length; A, B are 
the respective concentrations of acid and base forms of the reagent in the sample- 
reagent mixture; Cr is the total reagent concentration in the sample-reagent mixture; 
ej^ , £g are molar extinction coefficients at wavelength /I/ for A, B, respectively; and 
ODR^l is the optical density ratio as defined in Equation (5). 
[0005] The pH of a sample can be determined using spectral measurements by 
substituting Equation (6) in Equation (1). Because the acid and base form 
concentrations appear only as a ratio in Equation (1) and the absolute concentration of 
the reagent does not appear in Equations (1) or (5), the pH calculation is independent 
of the reagent concentration and the volume of the reagent added to the sample. If the 
ionic strength of the sample is at least an order of magnitude greater than the reagent 
concentration in the reagent-sample mixture, then the ionic strength of the sample and 
the activity coefficients (see Equations (7) and (8) below) are independent of reagent 
concentration. Accordingly, the only requirements for reagent addition are that the 
reagent concentration be (a) within a range where Beer's law is satisfied, (b) below an 
upper limit depending on the buffer strength of the sample beyond which the addition 
of the reagent could alter the sample pH, and (c) high enough to allow a good signal- 
to-noise ratio. 
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[0006] The molar extinction coefficients for the acid and base forms are obtained by 
calibration using solutions having pH values wherein the reagent exists completely in 
either its acid form or its base form. Alternatively, if the total reagent concentration is 
known very accurately, pH may be determined by measuring absorption at a single 
wavelength. However, because small errors in absolute concentration (Cr) can cause 
large errors in pH calculation, this method may not provide an accurate pH 
measurement. The pH of a sample may also be calculated using a continuous spectral 
scan in the relevant region and applying regression analysis to determine the base to 
acid ratio. 

[0007] From Equation (1), as pH values move away from the pKa values of a given 
reagent, the acid or the base fraction of the reagent becomes very small. Due to the 
low signal-to-noise ratios, the error in the pH measurement increases as pH values 
move away from the pKa value. Thus, for example, for a pH value 2 units lower than 
the pKa value, the fraction of the base form of the reagent is only 1%. Consequently, 
the OD corresponding to the base form peak wavelength is very low, resulting in 
inadequate pH accuracy. With single reagent indicators, the typical range of pH 
measurement is limited to about 1 to 1.5 units on either side of the reagent's pKa 
value. Beyond this range, a different reagent with a more appropriate pKa value must 
be used. Accordingly, the spectroscopic technique is less flexible to implement 
because knowledge of the sample's pH range is required a priori so that the 
appropriate reagent indicator is selected. 

[0008] A commonly used pH indicator uses a mixture of reagents to extend the range 
of pH measurement. Visual observation of color allows pH determination to within 1 
unit (see Vogel, "Text-Book of Quantitative Inorganic Analysis" 3"^"* Edition, Chapter 
1.30, page 59 (1961), incorporated by reference herein in its entirety). However, it is 
difficult to obtain highly accurate measurements using spectroscopic techniques 
because the visible spectrum of the mixed reagent is generally a cumulative addition 
of the spectra of individual reagents. Unless the individual spectra are well resolved, 
it is difficult to invert the fraction of each form for accurate pH calculation. 
[0009] Additional factors must be considered when performing pH analysis 
downhole. For example, only a limited number of reagents (pH indicators) can be 
transported downhole, the nature of the sample cannot be determined a priori, and 
only limited spectral analysis can be performed downhole. Accordingly, it would be 
useful to have a broad pH indicator that allows for simple, yet accurate, pH 
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determination. Further, it is difficult to control the amount of reagent added to the 
sample under investigation in the downhole environment. Accordingly, a pH 
measurement that is not dependent on reagent concentration would be preferred. 
[0010] Accordingly, it is an object of the present invention to provide an indicator 
mixture that retains the advantages of the single reagent method and that is effective 
over a broad range of pH values. 

[0011] It is a further object of the present invention to provide a spectroscopic 
technique to determine the pH of a sample with the accuracy levels comparable to 
single reagent spectroscopy over a broad range of pH values, 
[0012] It is yet another object of the present invention to provide a reagent mixture 
suitable for use in a downhole environment. 
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SUMMARY OF THE INVENTION 
[0013] The present invention relates to the use of a reagent mixture that allows pH 
measurement over a broader range of pH than available using conventional 
spectroscopic techniques. In one embodiment of the present invention, a fluids 
analyzer for determining the pH of formation fluid in a region of earth formation 
surrounding a borehole is disclosed comprised of: a reagent container having a 
mixture of two or more reagents wherein the reagent mixture is capable of detecting 
either: (1) a pH range broader than each reagent individually or (2) pH more 
accurately than each reagent individually; and means for mixing formation fluid with 
the reagent mixture downhole. The fluids analyzer may further include a spectral 
analyzer means, coupled to receive a mixture of formation fluid and the reagent- 
mixture from the mixing means for analyzing the optical density of the mixture of 
formation fluid and reagent mixture in two or more spectral channels. The reagent 
mixture is capable of measuring pH to an accuracy of about 0.5 units or better. 
Typical downhole spectral analyzers use spectral channels having bandwidths of 
approximately 10-20 nm. 

[0014] A second embodiment of the present invention is a method of making a 
reagent mixture for determining the pH of a sample comprising: (a) identifying a 
target pH range; (b) mixing known relative quantities of two or more reagents to 
create a reagent mixture, wherein the reagent mixture is capable of detecting either: (i) 
a pH range broader than each reagent individually or (ii) pH to a better accuracy than 
the reagents individually; (c) optimizing the reagent mixture to identify one or more 
sets of reagents to be mixed that satisfy the targeted pH range; and (d) characterizing 
at least one set of reagents. 

[0015] The sets of reagents may be optimized by identifying one or more sets of 
reagents by carefully selecting the thermodynamic acid dissociation constants of the 
reagents, adjusting the relative quantities of the reagents, or adjusting (the bandwidth 
^and/or wavelength of) the spectral channels. This optimization process may be 
performed by trial and error, or experimentation. It is noted that in some cases 
optimization may be generally viewed as an error minimization technique; however, 
optimization may be performed to improve other parameters of interest, such as 
maximizing sensitivity of pH measurement or signal-to-noise ratio. Alternatively, 
forward models and algorithms may be developed based on thermodynamic acid 
dissociation constants, relative quantities (or fractions), and the spectral channels. 
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Once "optimal" reagent sets are identified, at least one of the sets should be 
characterized to determine the detectable pH range of the mixture. Further, an error 
analysis should be performed to determine the accuracy of the mixture to measure pH. 
It is further noted that the absolute concentration of the reagent mixture may be 
optimized to improve signal-to-noise ratio and reduce the effects of error. 
[0016] In a third embodiment, a method for determining the pH of formation fluid in 
a region of earth formation surrounding a borehole is disclosed comprising: (a) storing 
a reagent mixture in a reagent container coupled to a fluids analyzer, wherein the 
reagent mixture is capable of detecting either: (i) a pH range broader than each 
reagent individually or (ii) pH to a better accuracy than each reagent individually; (b) 
positioning the fluids analyzer downhole; (c) drawing formation fluid into the fluids 
analyzer; (d) mixing formation fluid and the reagent mixture; and (e) analyzing the 
optical density of said mixture of formation fluid and reagent mixture in one or more 
channels to determine the pH of the formation fluid. If only one spectral channel is 
used, the absolute concentration of the reagent mixture should also be determined. 
However, if two spectral channels are used, then it is not necessary to determine the 
absolute concentration of the reagent mixture. 

[0017] Accordingly, the present invention provides a pH indicator mixture for 
spectroscopic techniques that is accurate over a broader pH range (as compared to 
conventional spectroscopic methods). As will be described in detail below, the 
method disclosed herein uses either (1) two or more spectral channels or (2) one 
spectral channel and accurate (absolute) reagent concentrations. If optical density 
(OD) measurements (or transmission, absorption or other spectral measurements 
which may be converted to OD) at two spectral channels are used, reagent 
concentration in the sample-reagent mixture is not needed; instead the relative 
quantities (fractions) of the reagent in the mixture may be used. Furthermore, the 
reagents need not have well resolved individual spectra. In fact, the accuracy of the 
method is enhanced if the spectra of the acid and base forms of one reagent are 
substantially similar to the spectra of the acid and base forms of the other, 
respectively. Because there is a vast pool of reagents to choose from, it is relatively 
simple to pick reagents that show similar direction of color change with pH though at 
different pKa values. Accordingly, the present invention includes algorithms to 
optimize for the choice of reagents (i.e., their pKa values) based on the pH 
measurement range of interest. Further, the individual reagent fractions in the mixture 
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may be selected to obtain an optimized mixture over the pH range of interest. It is 
noted that while the discussion that follows relates to the use of two reagents, the 
method may be easily extended to three or more reagents. 

[0018] For the purposes of the present invention, the term "optical techniques'* is 
intended to apply to a broad range of techniques, including but not limited to, 
spectroscopic, colorimetric, optical, etc. It is further noted that for the purposes of the 
present invention, the terms "pH sensitive reagent," "reagent," and "dye indicator" are 
used interchangeably. 

[0019] Further advantages and novel features of the invention will be set forth in the 
description which follows. 



7 



BRIEF DESCRIPTION OF THE DRAWINGS 
[0020] Figure 1 is a plot showing the spectra of a three-reagent mixture in buffer 
solution with pH varying from 3 to 9. 

[0021] Figures 2(a) and (b) are plots: (a) comparing model predicted pH from 
experimental ODR values of Figure 1 with true buffer pH values (445 and 570 nm 
wavelengths) and (b) showing the error in pH measurements as a function of pH. 
[0022] Figures 3(a) and (b) are plots: (a) showing a scenario where the acid form 
spectra of one reagent is identical to the base form spectra of the other and vice versa 
resulting in the pH being a non-unique function of ODR and (b) showing that pH is a 
unique function of ODR even when the acid and base forms of the two reagents have 
identical spectra. 

[0023] Figure 4 is a plot showing the error estimation for single and two-reagent 
mixtures, allowing estimation of target pH intervals for any desired accuracy level 
(see Table 2 below). 

[0024] Figure 5 is a plot showing the effect of changing the power p in the objective 
function [see Equation (30) below], 

[0025] Figure 6 is a plot showing that the optimization of the mole fractions in the 
three-reagent mixture yields better accuracy over the pH interval 4-9 as compared to 
using equal mole or weight quantities. 

[0026] Figure 7 is a plot showing the change in optimum reagent quantities with 
change in the target pH interval. 

[0027] Figure 8 is a plot showing the determination of the optimum pKa values for an 
equal weight mixture of three-reagents for various target pH intervals by optimizing 
objective function Equation (30) below. 

[0028] Figures 9(a) and (b) are schematics describing the overall optimization 
algorithm. 

[0029] Figure 10 is a schematic showing one application of the present invention. 
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DETAILED DESCRIPTION 
[0030] The present invention relates to a spectroscopic technique accurate over a 
broader pH range (as compared to conventional spectroscopic methods) using reagent 
mixtures. As will be described in detail below, the method disclosed herein uses 
either (1) two or more spectral channels or (2) one spectral channel and accurate 
reagent concentrations. If two spectral channels are analyzed, reagent concentration 
in the sample-reagent mixture is not needed; instead the relative quantities (fractions) 
of the reagent in the mixture may be used. Furthermore, the reagents need not have 
well resolved individual spectra. In fact, the accuracy of the method is enhanced if 
the spectra of the acid and base forms of one reagent are substantially similar to the 
spectra of the acid and base forms of the other, respectively. Because there is a vast 
pool of reagents to choose from, it is relatively simple to pick reagents that show 
similar direction of color change with pH though at different pKa, Accordingly, the 
present invention includes algorithms to optimize for the choice of reagents (i.e., their 
pKa values) based on the pH measurement range of interest. Further, the individual 
reagent fractions in the mixture may be selected to obtain an optimized mixture over 
the pH range of interest. It is noted that while the discussion that follows relates to 
the use of two reagents, the method may be easily extended to three or more reagents. 
[0031] Equations (1) through (6) may be rewritten for two reagents as follows: 



pH = pKa, +log^ + log4^ 


(7) 


pH = pKa^+ log + log ^ 


(8) 




(9) 




(10) 




(11) 




(12) 


OD^^ = eiflA^ + e^flB^ + e^lA^ + e^lB^ 


(13) 




(14) 
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where Ay, Bi are acid and base form concentrations, respectively, of reagent 1 in the 
sample-reagent mixture; A2, B2 are acid and base form concentrations, respectively, of 
reagent 2 in the sample-reagent mixture; C7/, C72 are total concentrations of reagent 1 
and reagent 2, respectively, in the sample-reagent mixture; /; is the mole fraction of 
reagent 1 in the reagent mixture; pKaj^ pKa2 are dissociation constants for reagents 1 
and 2, respectively. 

[0032] From Equations (7) through (14): 



ODRi^ = 



XI ( 



A2 



z 



-Bi 



B, 



IQ-(P»-Pf<ai) I ^Bi 
\ ^Ai J 



for/ =1,2 (15) 



where 



_ i-no-(^^'^^--) (1-/J 

5, \^\Q<p^-p^-^) f, ^^^^ 



Hence: 



ODRi^ =f(pH,t,pKa,,pKa^,f,) (17) 



where 

c-/o^l A2 _A1 _>12 _A2\ ^1 

^ '*^B1 »*^A1 '^fll *^A2^^B2^^A2*^B2 J V^O/ 



[0033] From Equation (17), which is an analytical equation, once the reagents for the 
two-reagent mixture are chosen, the relative quantities, and wavelengths (width of the 
spectral channel) are selected, the ODR becomes a function only of the pH. It is 
independent of absolute reagent concentrations or volume of reagent mixture added to 
sample, and therefore retains the advantage of the single reagent measurement. The 
mixture may be characterized in terms of ODR and pH such that plots and look-up 
tables may be generated . Alternatively, Equation (17) may be inverted to express pH 
as an analytical function of ODR (such as by using a program like Mathematica by 
Wolfram Research, Inc.): 
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pH = fiODR^^ . e, pKa, , pKa^ , /, ) 
[0034] Equations (7) through (19) can be extended to a system of n reagents: 



(19) 



A2 ( 



ODRi^ = 



1 V 









( 


£^'^ 










^Ai J 






^B\ 



fori=l,n (20) 



where: 



Hence: 



c-/c-^l ^>i2 ^^1 _yll ^/12 ^/ll -A2 ^A2 ^ 

"V'^v^l '^Sl ''^Al '*^i42'^B2'^/i2 »^fi2 ••**^>\n »^An / 



for/=l,n (21) 



(22) 



(23) 



(24) 



It is noted that Equations (22) and (24) are analytical expressions. 
[0035] Figure 1 shows experimentally measured spectra for a three-reagent mix in pH 
buffers 3-9. Arrow A indicates decrease in OD of acid forms of reagents as buffer pH 
increases from 3 to 9. Arrow B indicates increase in OD of base forms of the reagents 
as buffer pH increases from 3 to 9. Reagent properties of these mixtures are 
summarized in Table 1. 

[0036] Figure 2(a) compares model predicted pH values calculated using measured 
ODR values of Figure 1 with true pH values of the buffer solutions. Figure 2(b) plots 
the experimental error at various pH values. The reagent mixture of the example of 
Figure 2(b) is accurate over the range of buffers used (pH 3-9) with the errors within 
0.06 pH units. Accordingly, the reagent mixture of the present invention is accurate 
over a broader pH range than the traditional single reagent indicators, which typically 
work over 2-3 pH units (see, for example. Table 2 below). 
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Table 1: Summary of reagent properties used in simulations 





Phenol Red 
(reagent 1) 


Chlorophenol Red 
(reagent 2) 


Bromophenol Blue 
(reagent 3) 


Acid 


Base* 


Acid 


Base* 


Acid 


Base* 


s at 570 nm (A2) 


108 


37975 


58 


5A2A1 


378 


46859 


e at 445 nm (Xi) 


17916 


3352 


18136 


1985 


21711 


1981 


Mol. wt. 


354.38 


376.36 


423.28 


445.26 


669.98 


691.97 


pKa 


7.79 


6. 


11 


4. 


LI 



* The base form is a sodium salt of the reagent. 



[0037] The pH as calculated from Equation (24), however, is not necessarily always a 
unique function of ODR. Figure 3(a) is a plot showing two reagents (pKa]=7J9, 
pKa2=^6Al,fj=0,5) that are exact mirror images of each other; in this example, the 
spectrum of the acid form of one reagent is identical to the spectrum of the base form 
of the second reagent and vice versa (f^i = ^bi'^^ai ^ 7)- ^or each pH value 

there is a unique ODR. However, the inverse is not true; for a given ODR, there are 
two pH values and the method fails. By choosing reagents that show conrmion 
direction of color change with pH (but at different pKa values), the sensitivity of the 
pH to ODR may be increased to obtain a unique solution. 
[0038] Figure 3(b) shows the pH vs. ODR plot for two reagents (pKaj=7 J9, 
pKa2=6.11,f]=0,5) with identical spectra for its acid and base forms 
(^M = ^^^2 '^51 = ^B2 for J)' Here, the present invention works well. These results 
are contrary to conventional thinking where good spectral resolution would be 
preferred to calculate pH by inverting the spectra of the mixture for acid and base 
fractions of each reagent. 

[0039] Like the single reagent measurement, the mixture of the present invention is 
independent of total amount of reagent added to the sample. Accordingly, the present 
methodology is easy to implement in batch or flow injection mode for measuring pH 
of aqueous samples. The method uses (1) a minimum of two spectral channels if 
concentration independence is desired or (2) one spectral channel and accurate 
reagent concentrations. The latter is very sensitive to absolute reagent concentrations 
and, hence, is potentially less robust. Alternatively, pH may be calculated using 
regression analysis, if spectral absorption (OD) data at more than two wavelengths or 
continuous scans are available. 
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Error Propagation and Analysis in pH Measurement 

[0040] As mentioned above, a single reagent is sensitive only over a limited pH range 
around its pKa value. The extent of this range is determined by the accuracy required 
in the pH measurement and the spectroscopic noise. An error analysis should be 
performed to determine the measurable range as a function of these two parameters 
for single reagents and mixed reagents. An error analysis can also be performed to 
optimize for one or more of: the choice of reagents, relative quantities of the reagents, 
or spectral channels. 

[0041] The method of error propagation is used to determine the standard deviation in 
pH measurement as a function of spectroscopic noise (standard deviation in OD) 
(Bevington et aL, "Data Reduction and Error Analysis for the Physical Sciences/' 3^^ 
Edition, Chapter 3, pages 36-41, (2003), incorporated by reference herein in its 
entirety). 

[0042] For x=f(u,v), the standard deviation in x ((Tx) as a function of standard 
deviations in u (<TM)and v (cTv) is given as: 



+ <7 



dx 
dv 



(25) 



[0043] For a given reagent mixture (fixed pKai, s,^) and for a given set of 
wavelengths (A/, A2), the error propagation algorithm may be applied to Equations (20) 

and (24): 



X2 



^2 



OD 



X2 



OD" 



(26) 



(27) 



[0044] Combining Equations (26) and (27), standard deviation in pH (op//) 
measurement is obtained as a function of spectroscopic noise at the two wavelengths 
Ay andA2: 
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_ dpH ,2 f ZOD^ . { foDl^Y 

'^'''^ ^ aoD^,^^ ^^^^^ ""t^J 

(28) 



[0045] The derivative of pH with respect to ODR may be determined using the 
analytical expression of Equation (24) (which may be determined by inverting 
Equation (22) (such as by using packages such as Mathematica by Wolfram Research, 
Inc.). 

[0046] Alternatively, because the expression is very long, the derivative from 
Equation (22) may be numerically calculated provided ODR is a continuous function 
of pH and differentiable at every point in the interval of interest. 



= 1 (29) 
dODRj^ dODRj^ 

dpH 



[0047] Figure 4 compares the standard deviation in pH as a function of standard 
deviation in OD (where spectroscopic noise is assumed to be 0.01 OD units and Cj-is 
2x10'^ M) for single reagents and for equal weight fraction mixture of the two 
reagents, PR and CPR from Table 1. This figure allows the measurable pH range to 
be estimated for a desired accuracy as a function of spectroscopic noise. Accordingly, 
the two-reagent mixture of this example allows measurement over a range longer than 
the single reagents mixtures (see also Table 2). It is noted that the addition of a third 
reagent will allow further extension of the measurement range. 



Table 2: Extending range of pH measurements with reagent mixtures for 0.1 pH unit 
desired accuracy assuming spectroscopic noise of 0.01 OD units 



Reagent 


Mole Fraction (/i.A /j) 


Ct, [Ml 


pH range 


PR 


1.0, 0.0, 0.0 


2x10 = 


6.6-8.7 


CPR 


0.0, 1.0, 0.0 


2x10-= 


4.7-7.0 


BB 


0.0, 0.0, 1.0 


2x10-' 


2.8-5.1 


PR-CPR 


0.55, 0.45, 0.0 (equal weight fractions) 


2x10"= 


5.1-8.2 


PR-CPR 


0.55, 0.45, 0.0 (equal weight fractions) 


4x10' 


4.8-8.6 


PR-CRP-BB 


0.42, 0.35, 0.23 (equal weight fractions) 


2x10' 


3.6-7.9 


PR-CRP-BB 


0.46, 0.26, 0.28 (optimized for pH 3.5 to 8) 


2x10' 


3.5-8.1 


PR-CRP-BB 


0.42, 0.35, 0.23 (equal weight fractions) 


4x10"* 


3.2-8.5 
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[0048] From Equation (28), assuming equal noise at two wavelengths 
( ^oD^* ~ ^oD^2 ~ ^oD )» ^pH increases linearly with god- Note also that while the pH 
calculation is independent of absolute reagent concentrations, the error in the pH 
measurement depends on absolute reagent concentration through the absolute OD 
values that appear in Equation (28). For a given reagent mixture with a fixed relative 
quantity of the individual reagents, the error in pH, ap//, decreases linearly with 
increasing absolute reagent concentration. However, as explained above, if 
concentration of the reagents added is very high, then concentration could affect the 
pH of the sample and introduce additional errors. Furthermore, the concentration may 
go beyond the range where Beer's Law is valid. 

Optimization of Reagent Mixtures 

[0049] Using the error estimation algorithm described above, the reagent mole 
fraction in the mixture may be optimized to minimize the error over a desired pH 
range. A possible formulation for the objective function for desired accuracy of 0. 1 
pH units, for example, over the interval pHi to pH2, could include: 

Minimize F(/,,/2,.../J=f; m;,(^7^^^ -O.l)" forall a^^>0.l (30) 

where m is the number of discrete points into which the interval between pHi and pHa 
is divided and is the weight applied to the error contribution at each point. Both 
and p may be chosen to adjust the contribution of the error at each pH value in the 
desired range. This method provides the flexibility of having different levels of 
accuracy within a broader constraint of 0. 1 pH units error. For example, by assigning 
a higher weight to the points in the interval between 5 and 6, one can optimize for 0.1 
unit accuracy over the entire range of 4 to 9 with a higher level of accuracy between 5 
and 6. 

[0050] Figure 5 shows the effect of varying power p with w/ equal to one when 
optimizing over a pH range of 4-9 for a three-reagent mixture (reagent properties are 
summarized in Table 1). At pH equal to 9, the deviation from 0.1 units is very high. 
By using a power of 1, the contribution of error at 9 to the objective function Equation 
(30) is 0.27 units and is about 9 times higher than the contribution of the error at pH 5 
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(0.03 units). However, if the power p is 0.2, the contribution to the objective function 
of the error at pH 9 (0.77 units) is only 1.54 times higher as compared to the objective 
function error at pH 5 (0.5 units). In the second case with power p equal to 0.2, the 
optimization algorithm gives nearly equal importance to minimizing the error at pH 5 
and pH 9, while in the first case the pH 9 error term dominates and the optimization 
algorithm ignores the error at pH 5. This is evident in the plots for the two cases. 
[0051] While this is one form of the optimization function, one skilled in the art 
would recognize that other algorithms may be used to optimize reagent mixtures. 
[0052] Standard minimization routines (such as the non-gradient simplex or amoeba 
DUMPOL algorithm in EVTSL package from Visual Numerics, Inc.) can be used for 
the optimization of reagent mole fractions in the mixture. 

[0053] Figure 6 shows the results of a three-reagent mixture optimized for pH range 
of 4 to 9 and compares with cases when equal weight and equal mole fractions are 
used, [Spectroscopic noise is assumed to be 0.0 1 OD units and total reagent 
concentration is 2x10"^ M. For Equation (30), p is 0.5 and is 1.] 
[0054] Figure 7 shows how the optimization algorithm (see Equation (30)) adjusts 
mole fractions of the reagents in the mixture as the pH range of interest changes. As 
the pH range (pH 4-9) is shifted to a smaller interval of pH 6-9, the fraction of the 
reagent 1 (pKa=7J9) increases, while that of reagent 3 (pKa=4Al) falls to zero. 
[Spectroscopic noise is again assumed to be 0.01 OD units and total reagent 
concentration is 2x10"^ M. For Equation (30), p is 0.5 and is 1.] 
[0055] As discussed above, it is possible to obtain reagent combinations where the pH 
is a non-unique function of ODR, To prevent this, it is recommended that when a set 
of reagents are chosen, the mixture is forward modeled with known relative quantities 
or mole fractions of reagents to calculate pH as a function of ODR, The plot can be 
visually inspected or algorithms may be used to inspect the pH values to check for 
uniqueness. Generally, it is expected that when reagents in the mixture have different 
pKa values but their acid and base forms show similar color transformations with pH 
(for example, colors of all reagents changing from yellow to blue/red as they go from 
acid to base), the reagent combination should provide a unique pH vs. ODR plot. 
[0056] It is also possible to optimize the reagent selection process (with a final check 
being done as discussed above for uniqueness). Thus, for a fixed mole fraction, pKa 
values of the reagents in the mixture may be optimized so as to maximize accuracy in 
the pH interval of interest. Figure 8 lists optimum sets of pKa values of reagents in a 
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three-reagent mixture for different pH intervals of interest. [For Equation (30), p is 
0.5 and is 1.] For this example, the three-reagents had equal weight fractions, and 
epsilon values were fixed as provided in Table 1. 

[0057] Reagents may also be chosen to allow pH measurement over discontinuous pH 
intervals of interest. Thus, in a two-reagent mixture of phenol red (pKa=7J9), for 
example, and bromophenol blue (pKa=4Al) where the pKa values are separated by 
more than 3 units, the pH intervals on the low end (2.8-5.1 units) and the high end 
(6.6-8.7) may be simultaneously targeted without any sensitivity for regions in 
between. 

[0058] The general method of making and optimizing the pH indicator mixture is 
described in Figures 9(a) and (b). As shown in Figure 9(a), first a target pH range is 
selected 900. Then, known relative quantities of the two or more reagents are mixed 
910. The reagents are selected such that when mixed, the reagents cover a broader pH 
range or improved accuracy than if the reagents were used individually. It is noted 
that the present method does not require that absolute concentrations of the reagents 
be known; knowledge of relative (or fractional) quantities is sufficient. The reagent 
mixture is then optimized (for example, by trial and error, experimentation, or using 
forward modeling) to identify one or more sets of reagents that satisfy the target pH 
range 920. Once optimized reagent sets are identified, one or more of the reagent sets 
may be characterized to be used for spectroscopic pH measurement 930 (such as by 
developing functions of pH, OD, ODR, error deviations, etc.). 

[0059] A non-limiting preferred algorithm is shown in Figure 9(b). As in Figure 9(a), 
a pH target range is identified 900 and known quantities of two or more reagents are 
mixed 910. Boxes 920a-920f describe a preferred method of optimizing the reagent 
mixtures. First, a forward model is developed 920a based on the three variables: (1) 
the thermodynamic acid dissociation constants which are dependent on the reagents 
selected; (2) the relative quantities of the reagents (it may be simplest to begin with 
equal quantities of reagents); and (3) the identification of two or more spectral 
channels wherein the reagents show elevated levels of optical densities. It may be 
preferred to identify common channels wherein all of the reagents show elevated 
levels of optical density. Next, an algorithm is developed 920b based on the forward 
model of 920a to identify sets of reagents having optimized thermodynamic acid 
dissociation constants. (It is noted that reagents having the exact thermodynamic acid 
dissociation constants as optimized may not be available. Accordingly, an 
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"optimized" reagent set will include reagents have thermodynamic acid dissociation 
constants within acceptable range of the optimal value. Accordingly, numerous 
acceptable reagent sets may be identified.) 

[0060] The procedure is then repeated for the optimized set of reagents 920c-920d 
using the optimized thermodynamic acid dissociation constants, relative weight 
fractions, and spectral channels. However, now the reagent sets are optimized to 
determine the optimal relative weight fractions for the reagent sets identified in 920b. 
The procedure is repeated once more for the optimized set of reagents 920e-920f 
using the optimized thermodynamic acid dissociation constants, the optimized relative 
weight fractions, and spectral channels. But, now the reagent sets are optimized to 
determine the optimal spectral channels for the reagent sets identified in 920d. 
[0061] It is noted that each level of optimization need not be performed. One may 
chose to optimize for one or more (but not necessarily all) of: the thermodynamic 
equilibrium constants, relative quantities, or spectral channels at any time as shown by 
the dotted arrow to the left. At any point in the process, the optimized reagent sets 
may be checked for a unique correlation of pH as a function of ODR 930a and 
characterized for later use for spectroscopic pH measurement 930b. It is noted that 
the absolute concentration of the reagent mixture may also be adjusted to improve 
error (i.e., the signal to noise ratio.) 

Apparatus Example 

[0062] The present invention may be implemented in a laboratory or in various 
downhole fluid analysis tools. For example, the apparatus described in commonly 
owned co-pending United States Patent Application Serial No. 10/667,639 filed 
September 22, 2003, entitled "Determining Fluid Chemistry of Formation Fluid by 
Downhole Reagent Injection Spectral Analysis" (incorporated by reference herein in 
its entirety) is a preferred implementation of the present reagent mixture. 
[0063] One non-limiting embodiment of the present invention, as shown in Figure 10, 
is a wireline formation tester 1015, including fluids analyzer 1020. The formation 
tester is shown downhole within borehole 1012 of formation 1010 suspended by 
logging cable 1017. Logging cable 1017 also couples the formation tester to surface 
system 1018. Fluids analyzer 1020 is located within formation tester 1015. Means 
for extracting fluid from the formation, in this case a probe 1021, is shown extending 
from the formation tester 1015. Major elements of fluids analyzer 1020 are fluid 
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flow-line 1030, reagent container 1031, and spectral analyzer 1041. This embodiment 
further uses a flow-line pump (dual chamber piston pump) 1051 and main pump-out 
line 1061. It is noted that this figure is intended to depict a non-limiting embodiment 
useful for deploying the present invention. Other suitable elements may be included 
as dependent upon the specific application. For example, other configurations may be 
used to extract fluids from the formation, such as in a producing well or a monitoring 
well where a pump is used to obtain fluid samples. Accordingly, while Figure 10 
depicts a wireline application, the present invention may be deployed in other 
scenarios, such as permanent installations or a drilling environment. 
[0064] While the invention has been described herein with reference to certain 
examples and embodiments, it will be evident that various modifications and changes 
may be made to the embodiments described above without departing from the scope 
and spirit of the invention as set forth in the claims 
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